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The paper deals with the problem of integration of equations of motion in nonholonomic systems. By means of well-known 
theory of the differential equations with an invariant measure the new integrable systems are discovered. Among them 
there are the generalization of Chaplygin's problem of rolling nonsymmetric ball in the plane and the Suslov problem 
of rotation of rigid body with a fixed point. The structure of dynamics of systems on the invariant manifold in the 
integrable problems is shown. Some new ideas in the theory of integration of the equations in nonholonomic mechanics 
are suggested. The first of them consists in using known integrals as the constraints. The second is the use of resolvable 
groups of symmetries in nonholonomic systems. The existence conditions of invariant measure with analytical density 
for the differential equations of nonholonomic mechanics is given. 



1. Introduction 

The integration theory of equations of motion for mechanical systems with nonholonomic constraints 
isn't so complete as in the case of systems with holonomic constraints. This fact has many reasons. 
First, the equations of nonholonomic mechanics have a more complex structure than the Lagrange 
equations, which describe the dynamics of systems with integrable constraints. For example, a non- 
holonomic system can't be characterized by the only function of its state and time (cf. [1], ch. XXIV). 
Second, the equations of nonholonomic mechanics have no invariant measure in the general case (a sim- 
ple example is given in section 5). The point is that nonholonomic constraints may be realized by 
action of complementary forces of viscous anisotropic friction with a large viscosity coefficient ([3]). 
The absence of an invariant measure is a characteristic property of systems with friction. In limit, 
the anisotropic friction is compatible with the conservation of total energy. But asymptotically stable 
equilibriums or limit cycles may arise on the manifolds of energy levels (cf. [4]), and this is the reason 
for nonexistence of additional "regular" integrals of motion. 

The most popular method to integrate the equations of nonholonomic dynamics is based on the 
use of the available first integrals or the "conservation laws" : if a Lie group that acts on a position 
space preserves the Lagrangian and if vector fields that generate this group are the fields of possible 
velocities then the equations of motion have the first "vector" integral or the generalized integral of 
kinetic moment [5,6]. A number of problems of nonholonomic dynamics was solved by this method, 
among them, we note especially Chaplygin's problem on an asymmetrical ball rolling over a horizontal 
plane [5]. 

Attempts to generalize the Hamilton - Jacobi method to the systems with nonholonomic con- 
straints were non-effective as well as attempts to present the equations of nonholonomic dynamics in 
the form of Hamiltonian canonical equations. It turned out that with the help of the Hamilton- Jacobi 
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generalized method it is possible to find at most only some special solutions of the motion equations. 
This paper contains the detailed analysis of these questions. 

Another general approach to the integration of nonholonomic equations is based on the theory 
of Chaplygin's reducing factor ([5]): one try to obtain a change of time (different along different 
trajectories), such that the equations of motion are presented as Lagrange or Hamilton equations. 
Though such change exists in exceptional cases only, it allows to solve a number of new problems of 
nonholonomic dynamics (cf . [5] ) . Let us note that the equations of motion sometimes may be reduced 
to the Hamiltonian form by other reasoning (see section 5). 

The list of exactly solvable problems of nonholonomic dynamics isn't long: almost complete 
information may be found in the books [1,5,8]. In this work we present some new integrable problems, 
consider the characteristic features of behavior of nonholonomic systems' trajectories in the phase 
space, and propose some general theoretical reasonings on methods of integrating the equations of 
nonholonomic dynamics. 



2. Differential equations with an integrable measure 

Let us consider a differential equation 

x = f{x), x G R n (2.1) 

and let g t be its phase flow. Suppose (|2.1|) has an integral invariant with some smooth density M(x), 
i.e. for any measurable domain D C R n the following equation holds for all t 

M(x) dx = J M(x) dx. (2.2) 

g*(D) D 

Let us recall the well-known Liouville assertion: a smooth function / : R n — > R is a density of an 
invariant 

/*(.)* 

if and only if div(M/) = 0. If M(x) > for all x then Q2.2JI defines a measure in R n invariant 
with respect to the action of /. The existence of an invariant measure simplifies the integration of a 
differential equation; for example, in the case of n = 2 the equation is always integrable in quadratures. 
According to Euler, M is also referred to as an integrating factor. 



Theorem 1. Suppose system (|2.1jl with invariant measure (|2.2jl has n — 2 first integrals 
Fi, F n _2- Suppose the functions F\, F n _2 are independent on the invariant set E c = 
= {x E R n : F s {x) = c s , 1 s$ s < n- 2}. 

Then 

1) the solutions of (j2.1j) that belong to E c may be found by quadratures. If L c is a connected 
compact component of the level set E c and f ^ on L c then 

2) L c is a smooth surface diffeomorphic to a two-dimensional torus, 

3) it is possible to choose angle variables x, y mod 2ir on L c so that, after the change of variables, 
system (|2.1j) on L c would have the following form 

11 (2-3) 



${x, y y $(z, yy 

where A, /i = const, |A| + ^ and & is a smooth positive function 2n -periodical with respect to x 
and y. 
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Let us mention the main points of the proof. Since the vector field / is tangent to E c , differential 
equation l|2.1|) is bounded on E c . This equation on E c has an integral invariant 



/ 



Mda 

K-2 ' 



where da is the element of area of E c considered as a surface embedded into IR" , V n -2 is the (n — 2)-di- 
mentional volume of the parallelepiped in M n , the gradients of F\, . . . , F n -2 being its sides. Now the 
integrability by quadratures on E c follows from Euler's remark. The first conclusion of theorem 1 
(which was firstly mentioned by Jacobi) is proved by this reasoning. The second conclusion is the 
well-known topological fact: any connected, compact, orientable, two-dimensional manifold that ad- 
mits a tangent field without singular points is diffeomorphic to a two-dimensional torus. The third 
conclusion is, in fact, the Kolmogorov theorem on reduction of differential equations on a torus with 
a smooth invariant measure [9]. 

Equations (|2.3j) have invariant measure 

3>(x, y) \dx A dy. 

By averaging the right-hand sides of Q2.3|) with respect to this measure we get the differential equations 

u=£, v = £; „=-L J J $dxdy. (2.4) 



Proposition 1. Let <I>: T 2 — > R be a smooth (analytical) function. Then for almost all pairs 
(A, /i) € M there exists a smooth (analytical) change of angle variables x, y — > u, v that reduces (|2.3j) 
to (E2J). 

The proof is presented in [9,10]. Note that if ()2.3|) can(not) be reduced to (|2.4j) for a pair (A, /i) 
then the same is true for all pairs (xA, >c/j,), x 7^ 0. So, the property of reducibility depends on 

arithmetical properties of ^ that is referred to as the rotation number of the tangent vector field on 

T 2 = {x, y mod 2vr}. 

Proposition 2. Lei $(2;, y) = ^ V'm.n expi(mx + ny), </? min = Tp_ m _ n . If (|2.3j) may be reduced 
to (|2.4|) by a differ entiable change of angle variables u = u(x, y), v = v(x, y) then 

2 < 00. (2.5) 



E 



If the ratio — is rational then the torus T 2 is stratified into a family of closed trajectories. In this 

case the reducibility condition is equivalent to the equality of periods of rotation for different closed 
trajectories. 

In the general case (the Fourier decomposition of <3? contains harmonics) the points (A, n) £ ffi n 
with rationally independent (A, /i), for which series (|2.5j) diverges, are everywhere dense in R n . The 
questions of reducibility of (|2.3|) are discussed in [9]. The general properties of solutions of H2.3|) see 
in [10]. 



3. S. A. Chaplygin's problem 

Let us consider as an example the problem of rolling of a balanced, dynamically non-symmetric ball 
on a horizontal rough plane (see [5]). The motion of the ball is described by the following system of 
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equations in R 6 = R 3 {u} x IR 3 {7}: 

k + ojxk = 0, 7 + wx7 = 0; 

(3.1) 

k = Iuj + ma 7 x (w x 7). 

Let to be the vector of the angular rotation velocity of the ball, 7 the unit vertical vector, / the 
tensor of inertia of the ball with respect to its center, m the mass of the ball, and a its radius. These 
equations have the invariant measure with density 

M= 1 E = ||<%ll- ( 3 - 2 ) 

/ma 2 )" 1 - (7, (I + ma 2 E)- 1 -f) 

Taking into account the existence of four independent integrals F\ = (k, u), F2 = (k, 7), F3 = 
= (7, 7), i<4 = (k, k), we see that (|3.1j) is integrated by quadratures. Note that system of equa- 
tions ()3.1j) has no equilibriums on the non-critical level sets E c . Indeed, if 7 = 0, then u and 7 are 
linearly dependent. This fact implies the linear dependence of dF± and dF<i- The simplest case of 
integration by quadratures of equations (|3.1|) is the case of zero value of the constant in the "area" 
integral i*2- In elliptic coordinates £, rj on the Poisson sphere (7, 7) = 1 the equations of motion on 
the level E c are reduced to the following form 



Coefficients of the polynomial -P5 of the fifth order and the constant a depend on parameters of 
the problem and on the constants of the first integrals (see [5] for details). The variables £, rj ranges 
over different closed intervals a\ ^ £ ^ 02, b\ ^ r, ^ &2 where P5 is nonnegative. The uniformizing 
substitution 

^ I zdz yi _ 1 f zdz 



/TO *j /to 

/ 



(3.3) 



zdz -1 1 zdz 

y = t i , — v 



b 



/TO' V /TO 



introduces the angle variables x, y mod 2tt on E c , and the equations of motion take form (|2.3|) 



V- " 



(r 1 (x) - r," 1 (V)) ^(o -£(*)) (a -7,(2,)). 



Here £(x) and r,(y) are 27r-periodic functions of x and y arising as the inversions of Abelian inte- 
grals ipHj) . 

These equations imply 

Proposition 3. The rotation number of a tangent vector field on two-dimensional invariant 
tori in Chaplygin's problem is equal to the ratio of real periods of the Abelian integral 

z dz 



/TO 
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Remark. This assertion is true for integrable problems of dynamics of a heavy body with a fixed point 
defined by the system of Euler-Poisson equations (see [10]). Since the Euler-Poisson equations arc Hamiltonian, 
by the Liouville theorem, in integrable cases they always can be reduced to form (|2.4|) on two-dimensional 
invariant tori. It seems that equations l|3.4|) have no such property; inequality (|2.5[) is not fulfilled on all 
invariant non-resonant tori. 

Let us make a change of time t — > r by the formula 

dt = Oi a ~v) dr- (3.5) 

Equations (|3.4j) preserve their form but the variables x, y in function <3? are separated 

<S> = C 1 (x)-ir 1 (y). 
Proposition 4. Suppose that /i and A in (|2.3j) are nonzero and 

$ = &( X ) + $"(y). 

Then equations (|2.3j) can be reduced to form (|2,4|) by an invertible change of angle variables on T 2 . 
The proof is presented in [10]. Note that if <E> = &'(x) + $"(y), then series ()2.5|) 

v" 



' 2 



fjn 



+ 



nfjL 



2 



converges for all A, fi ^ 0. 

So, taking into account change of time Q3.5|) one may reduce (|3,4j) to the form 

^ = U, ^ = V, (3.6) 
ar ar 

where J7 and V depend on the constants of the first integrals only, and U, V ^ 0. This result can cause 
the temptation to use Chaplygin's reducing multiplier theorem: if we can reduce (|3.1|) using change of 
time ()3.5|) to the Euler-Lagrange equations of some variational problem (they are written, as well as 
the classical Euler-Poisson equations, in redundant variables) then according to the Liouville theorem, 
the equations of motion on the two-dimensional invariant tori in some angle variables u, v mod 2tt 
have form (|3.6|) . It is possible to show that this method does not lead to the goal. In conclusion, note 
that S. A. Chaplygin himself never considered the problem of the ball's rolling in connection with the 
reducing multiplier theory. 



4. A generalization of S. A. Chaplygin's problem 

We are going to show that the problem of rolling of a balanced, dynamically non-symmetric ball on 
a rough plane is still integrable (in the sense of section 1), if the particles of the ball are attracted 
by this plane proportionally to the distance. Since the center of mass of the ball coincides with its 
geometrical center, we can calculate the potential by the formula 

V(l) = | J{r, 7 ) 2 ^ = |(J 7 , 7 ), (4.1) 

where 7 is the unit vertical vector, r is the radius vector of particles of the ball, I is the tensor of 
inertia of the ball with respect to its center. The attraction forces generate the rotational moment 

— J r x (e(r, 7)7) dm = — e j (r, ^)(r x 7) dm = 7 x V' 1 = £7 x I7. 
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In order to obtain the moment of forces with respect to the contact point, it is necessary to add the 
moment of the combined force 

e j (r, 7)7 dm = el^j r dm, 7^7, 

which is equal to zero, since the center of mass of the ball coincides with its geometrical center. 

According to the theorem about the change of kinetic moment behavior with respect to the 
contact point (see [5] , [6] ) , the equations of rolling of the ball can be presented in the following form 

k + ujxk = e r yx Ij, 7 + uj x 7 = 0. (4.2) 

Theorem 2. Differential equations (|4.2|) are integrable by quadratures. 
Indeed, they have four independent integrals 

F 1 = (k, co) + e(J 7 , 7), F 2 = (k, 7), F 3 = (7, 7) = 1, 
F 4 = (k, k) - {A>y, 7), 

where elements Ai of a diagonal matrix A are expressed through the principal moments of inertia Ij 
by the formulae 

A\ = e(h + ma 2 )(h + ma 2 ), . . . 

Since equations (|4,2|) have the invariant measure with density ()3.2j) . they are integrable by theorem 1. 
It would be interesting to integrate this equation explicitly and test if proposition 3 remains true for 
equations (|4,2j) . 

Note that the problem of rotation of a body about a fixed point in an axisymmetric force field 
with potential (|4.1j) is also integrable ([1]). In addition to the classical integrals F±, F2, F3, there is 
the integral -F4, where one must put a = 0. This integral was found independently by Clebsh in the 
problem on motion of a body in an ideal fluid and by Tisseran, who investigated rotational motion of 
heavenly bodies. 

5. G. K. Suslov's problem and its generalization 

Following G. K. Suslov ([11], ch. 53), we consider the problem of rotation about a fixed point of a body 
with the nonintegrable constraint (a, u) = 0, where a is a vector that is constant in the moving frame 
of reference. Suppose that the body rotates in an axisymmetric force field with the potential V{^). 
Following the method of Lagrange multipliers, we write down the equations of motion ([11], ch. 46): 

Iuj + uj x Iuj = 7 x V' + Xa, 7 + ^x7 = 0, (a, u) = 0. (5.1) 

Using the constraint equation (a, uj) =0, the Lagrange factor can be expressed as the function 
of uj and 7 

A = -(a, r\luj xuj) + r 1 ^ x V!j))/{a, r x a). 
Equations (|5.1j) always have three independent integrals: 

F 1 = (Iuj, uj)/2 + V(j), F 2 = (7, 7), F 3 = (a, uj). 

For real motions, F 2 = 1, F 3 = 0. In this case, we can reduce the problem of integration of 
equations (|5.1|) to the problem of existence of an invariant measure (the existence isn't evident) and 
the fourth independent integral. 

Proposition 5. // a is an eigenvector of operator I, then the phase flow of system (j5.1|l 
preserves the "standard" measure in M 6 = R 3 {uj} x E 3 {7}. 

To prove the proposition we have to verify the following fact: the divergence of the right-hand 
side of (|5.1|) is equal to zero as la = fia. 
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G. K. Suslov has considered a particular case of the problem, when the 
body is not under action of exterior forces: V = 0. In this case the first 
equation of (|5.1|) is closed relatively to lj. We can show that it is integrable by 
quadratures (see [11], ch. 53). The analysis of these quadratures shows that if a 
isn't an eigenvector of the inertia operator, then all trajectories u(t) approach 
asymptotically as t — > ±oo to some fixed straight line on the plane (a, uj) = 
= (see Fig. 1). Consequently, the equation with respect to uj and complete 
system ()5.1|) have no invariant measure with continuous density. In this case 
theorem 1 isn't applicable, so, the question about the possibility to find the 
vector j(t) by quadratures remains open. But if la = fia then equations ()5.1I) 
have the additional integral: the value of the kinetic moment is preserved 




Fig. 1 



F A = (Iuj, Iuj). 

Equations (|5.1|) are integrable by theorem 1. However, this possibility may be easily realized directly. 
It seems that in the most general case, the existence of an invariant measure is connected with the 
hypothesis of proposition 5: la = \xa. From now on, we suppose that this equality is fulfilled. 

Now suppose that the body rotates in the homogeneous force field V = (b, 7). If (a, b) = 0, 
then equations (|5.1|) have the integral 

F 4 = (Iu, b) 

consequently, they are integrable by quadratures. This case was indicated by E. I. Kharlamova in 
her work [12]. We are going to consider an "opposite" case, when b = ea, e / 0. Without loss of 
generality we can assume that the vector a has the components (0, 0, 1). Taking into account the 
equation U3 = 0, we obtain that two first equations 1)5.1)) have the following form 

hdJi = £72, = -£71; w = (^1) ^2, ^3)- 



Therefore I\U)i = £72, = — £7i- 

Using the Poisson equations 71 = — W273, 73 = ^173 we get 

I\UJi = £73(^1, I2U2 = SJ2^2- 



(5.2) 



The energy integral 



(hcof + / 2 w|)/2 + £73 = h 



makes it possible to express 73 through lo\ and u>2- After that, equations 1)5. 2 Jl may be rewritten as 
the Lagrange equations 

h^i - 75— ^ -77157- - 7T7 l*- 1 ' 2 ); 



duJi dt duJi duji 



l\Co\ + 1^(02 



L = T V , 



V=±lh 



If IiUjj + I 2 ooi\ 



These equations have the energy integral T + V. For real motions its value is equal to £ 2 /2. Let us 
emphasize that unlike the reducing multiplier theory our reduction of equation 1)5.1)1 to Lagrange (or 
Hamilton) equations doesn't require the change of time (cf. [11]). 

The change ijWj = ki corresponding to the transition from the angular velocity to the kinetic 
moment reduces the considered problem of rotation of a body to the problem of motion of a material 
point in the potential force field 



kj 



-w (z-i 2) v-Uh fc^ 1 + fcfo- y 



(5.3) 
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At I\ = I2 we have the motion of point in a central field. This motion corresponds to the well-known 
integrable "Lagrange case" of the generalized Suslov problem. As well as in Lagrange's classical prob- 
lem of a heavy symmetric top, the equations of motion are integrable in this case in elliptical functions 
of time. If Ji 7^/2) then the equations apparently have no additional analytical integral independent of 
the energy integral. The following observation confirms this supposition. Put formally I\ = —I 2 = 1. 
Then at h = equations (|5.3f) practically coincide with the equations of the Young- Mills homogeneous 
two-component model, non-integrability of which is established in [14]. 

If the value of h is fixed, the point moves in the area defined by the inequality V ^ e 2 /2. For 
different h, these areas are shown in Fig 2. The trajectories of vibrational motions, when one of the 
components k\ or ki becomes zero, are especially interesting. These motions are expressed through 
elliptical functions of time. 

One more case of integrability of Q5.1|) is given by 

Theorem 3. Suppose la and the potential V(7) is 
defined by formula Ij4.1j) . Then equations (|5.1|) are integrable 
by quadratures. 





h > e 



Fig. 2 



Proof. 

f l< e Let us show that equations (|5. 1|) have the Clebsch-Tis- 

seran integral 



\{Ai, 7), 



A = el' 1 det /. 



Indeed, 

F A = (Iu, 7 x e/7) + \{a, Iu) + (^7, u x 7) = 

= (to, 1(7 x e/7)) + Xfi(a, uj) + (to, 7 x A'y) = 

= (uj, 1(7 x e/7) + 7 x A'y) = 0, 

since ^(7 x J7) = — (7 x i" _1 7)det7. To complete the proof, we have to take into account the 
conclusion of proposition 5 and to use theorem 1. 

Let us show, how one can explicitly integrate equation ()5.1j) . For definiteness, let a = (0, 0, 1) 
and e > 0, -Z3 ^ Ji, I3 ^ I2. Then (|5.1|) may be presented as the following closed system of four 
differential equations: 

hibi = e(ii - I 2 )7273, ^2^2 = e(7 3 - 11)7173, 



71 



-^273, 72 =0^73, 7! = 1 — 71 — 7I . 



Let us introduce the new time r by formula dr = 73 dt and denote the differentiating with respect 
to r by prime. Then equations of motion take form of a linear system with constant coefficients 

Iiu/i = e(J 2 - 13)72) 7'2 = w i> 
They can be presented in the equivalent form 



Y'l + A?7i = 0, 



Put 



A? 



arctg 



A171 

U>2 ' 



7" 2 + A^72 = 0, 

\l = {h-h)/h- 

A272 



ip 2 = arctg ■ 
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These variables are angle variables on two-dimensional invariant tori with 
Consequently, 

<j>i = \i/$, ip2 = A 2 /$; $ = (1 - c\ sin 2 c^i - c| sin 2 ip 2 ) 



-1/2 



The constants c\ and c 2 (c 2 + c 2 . ^ 1) can be expressed as functions of constant values of the energy 
integral and Clebsch - Tisseran integral. The remarkable property of this problem is the fact that 
the ratio of frequencies A1/A2 is independent of initial data and depends only on the constants of 
parameters of the problem. Consequently, if the number 



< {h-h)h 



is rational, then all solutions are periodic; otherwise practically all trajectories aren't closed (except 
degenerated motions, when 71 = or 72 = 0). Let y s (0) = a s . Then 

r 

dx 



1 — c 2 sin 2 (Aix + ai) — c 2 . sin 2 (A22; + a 2 ) 

If c\ = (or c 2 = 0) then ji and 72 (and consequently, U\, u>2, 73) are elliptical functions of time. 
This conclusion is true in the case Ai = A2 (i. e., when I\ = I2 or I3 = I± + I2) for all c\, 02- In the 
most general case the analytical character of the solutions is essentially more complex. In conclusion, 
note that series (|2.5|) diverges in this problem if the irrational ratio A1/A2 is approximated by rational 
numbers anomalously fast. 

Remark. Equations (|5.1(l are also integrable for potentials of the general form 

V{l) = ^(Cll7l + c 2272 + c 3373 + 2ci27l7 2 ). 

Using the change of time dr = 73 dt, the equations of motion are reduced to the linear system 

t V dV j 11 dV . f> T/ i 

7271 =~dTi ll2 =~W2' v - v k=i-^-yr 

In the general case, potential V does not have a simple physical interpretation. 

6. The first integrals used as constraints 

Let L(x, x, t) be a Lagrangian of a nonholonomic system that satisfies the following "regularity" 
condition: the quadratic form 

is positively definite. In particular, det 7^ 0. The constraints (not necessarily linear) are given 

by the equations 

fi(x, x,t) = ... = f m (x, x,t) = (6.1) 

with independent co- vectors 

dfi dfrn 
dx ' " " ' dx 
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The equations of motion can be presented as the Lagrange equations with the multipliers 

!(i)-i=£*<i. a-.-/.-* (6.2) 

8=1 

Proposition 6. If the functions L, /i, . . . , f m satisfy the above conditions, then a unique 
solution of (|6.2j) corresponds to every initial state that is permissible by constraints (|6.1|) . 

Indeed, under these suppositions the multipliers Ai, . . . , A m are smooth functions of x, x, t by 
the explicit function theorem. 

Now suppose equations (|6.2j) have a first integral F(x, x, t). We get the following 

Proposition 7. If the co-vectors ^Jr-, .... -rrr- are independent then x(t) is a solution 

ox ox ox 

of 1)6.2(1 with the constant value of the integral F = c if and only if this function is a motion of 
mechanical system with the Lagrangian L and constraints fx = . . . = f m = f m +i = c, where f m +\ = 
= F-c. 

The sufficiency is obvious for x(t) satisfies (|6.2[1 if we put A m +i = 0. On the contrary, let x(t) be 
a solution of a system of form (|6.2[) . where s ranges from 1 to m + 1. Let y(t) be the unique motion 
of system (16.211 with the initial data y(0) = x(0), y(0) = x(0). Evidently, F\ ,,,= c. The function y(t) 
as well as x(t) satisfies the equations of motion of the extended system with A m +i = 0. To complete 
the proof, it remains to use the conclusion of proposition 6. 

Let us discuss one of possible applications of proposition 7. Suppose fi are linear with respect to 
velocities and constraints <j6.1|) are non-integrable. If equations of motion have the linear integral F, 
then equations 

fl = ■ ■ ■ = fm = fm+l = (fm+l = F — c) 

may turn out to be completely integrable. In this case the study of motions that belong to the level 
set F = c is reduced to the investigation of some holonomic system. We do not have to integrate here 
the constraint equations, since the variables can be considered as redundant one's, and the equations 
of motion may be written as Hamilton equations in redundant variables (see [11], [15]). 

Let us consider as an example Suslov's problem in a homogeneous force field in the Kharlamova 
integrable case. The equations (a, uS) = and (Iuj, b) = form an integrable field of directions on 
the manifold of the rigid body positions (on the group SO(3)). Thus, Suslov's problem is reduced in 
this case to a system with one degree of freedom. Though the one-dimensional manifold of states of 
such system isn't closed in SO (3) in the general case. 

If the constraints are non-linear with respect to velocities, it is natural to use the energy integral 

H(x, x) = ^-x — L 
ox 

as the first integral 1 . 

For example, let us consider the Appell-Hamel system with the Lagrangian 

L = ^(x 2 + y 2 + i 2 ) + gz, g = const 

and the non-linear constraint 

x 2 + y 2 = k 2 z 2 , k = const / (6.3) 
(see [16] and [17]). By the energy integral 

hx 2 + y 2 + z 2 ) - gz = h 



1 The equation of motion has the integral of energy, if the constraints are homogeneous with respect to velocities and 
the Lagrangian does not depend on time explicitly. 
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and H6,3|) we get the equation of the "integrable" constraint 

z(l + k 2 ) , 
2 ~gz = h. (6.4) 

Consequently, the coordinate z changes with the constant acceleration g/(l + k 2 ). Excluding non- 
linear integrable constraint Q6.4|) (i. e., considering z as a known function of time) we get a more simple 
system with two degrees of freedom, the Lagrangian 

L = \(x 2 + y 2 ) 

and the constraint x 2 + y 2 = f(t), where / = k 2 z 2 is a known quadratic function of time. The further 
integration may be easily fulfilled. 



7. Symmetries of nonholonomic systems 

We suppose that the vector field v(x) ^ is a symmetry field of a nonholonomic system with La- 
grangian L(x, x) and constraints 

a;) = . . . = f m (x, x) = 0, 
if the phase flow of the differential equation 

dx 



dt V{X) 



preserves L and f\ , . . . , /„ 



Proposition 8. ^4 phase flow of a symmetry field converts solutions of a nonholonomic system 
to solutions of the same system. 

Proof. 

By the theorem on rectification of trajectories, the phase flow gf, in some local coordi- 
nates x±, . . . , x n is the following one-parameter group 

xi -> xi + s; x 2 -> » 2 , • • • ! x n —> Xn- 

With respect to these variables, L and fi do not depend on x±, consequently, the equations of motion 
do not contain this variable, too. This fact implies proposition 8. 

In the case of integrable constraints, the symmetry field corresponds to a linear with respect to 
velocities first integral of the equations of motion. It is not so in the case of nonholonomic systems. 

Proposition 9. If g% preserves the Lagrangian and v is the field of possible velocities, i. e. 

dx dx 
then the equations of motion have the first integral 

= const. 

ox 

This assertion ( "the Noether theorem") is discussed in [6], for example. 
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Theorem 4. Suppose the equations of motion ([6.2(1 have n — m first integrals f m +i, ■ ■ ■ > fn- If 

1) at points of the set E c = {/i = . . . = f m = 0, / m+ i = c m+1 , ...,/„ = £%} i/ie Jacobian 

d(h,...,f n ) 
dx\ , . . . , x n 

is nonzero, 

2) i/iere ercisi fields v±, i> n _i i/icrf are linearly independent at all points Ec and generate a 
solvable Lie algebra with respect to the commutation operation, while their phase flows gf,. preserve L 
and fi, f n , 

3) there are no vectors x = ^s v s( x ), A s G IR among solutions of the system of equations 

fl = ■ ■ ■ = fm = 0, fm+l = C m +1, . . . , fn = C n , (7.1) 

then solutions of (|6.2j) that belong to E c are found by quadratures. 

Remark. In some cases the existence of first integrals of nonholonomic systems can be established by the 
following observation. Let F(x, x) be the first integral of a "free" holonomic system with Lagrangian L. This 
function is an integral of a nonholonomic system with the same Lagrangian L and constraints /i = • • ■ = f m = 
in the case of 



fd 2 L\- 1 of s OF n , ^ ^ 



if fx = . . ■ = f m = 0. This condition of invariancy is fulfilled for the Clcbsch-Tisseran integral in Suslov's 
problem (theorem 3). Besides, it is fulfilled for the energy integral in the case of homogeneous constraints and 
dL 

for the Noether integral -jp- ■ v, if the field v is the field of possible velocities (proposition 9). 

Proof of theorem 4. 

By the explicit function theorem we obtain from ([7.1)1 that 

x = v c (x). (7.2) 

By conditions 2 and 3, the vectors v c , Vi, . . . , v n -\ are linearly independent at all points x. The phase 
flows g% convert solutions of (|7.2[) to solutions of the same equation (proposition 8). To complete the 
proof, it remains to apply the well-known Lie theorem on integrability by quadratures of systems of 
differential equations (see, for example, [18]). 

Let us consider as an illustrating example the problem of sliding of a balanced skate on horizontal 
ice. One can choose units of length, time and mass so that the Lagrangian would take the following 
form: 

L = \{x 2 + y 2 + z 2 ). (7.3) 

Here x, y are the coordinates of the point of contact, z is the angle of rotation of the skate. The 
constraint equation is 

f = x sin z — y cos z = 0. (7.4) 
The equations of motion have two first integrals 

x 2 + y 2 + z 2 = h, z = c. (7.5) 

The second one is obtained by using proposition 9 with the help of the vector field V3 = (0, 0, 1). 
By (|7.4[) and (|7.5jl we obtain the field Vh, c = (Vh — c 2 cos z, \Jh — c 2 sin z, c) . The fields vi = (1, 0, 0) 
and V2 = (0, 1, 0) are commuting symmetry fields. If c 7^ then vectors v^c, v\ and vi are linearly 
independent, consequently, in this case we can apply theorem 4. Let us emphasize that v\ and V2 do 
not generate the conservation laws. 
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Theorem 4 impose strict restrictions on the nonholonomic system. These restrictions can be 
weakened if we replace condition 2 by the condition 

2) for the fixed c = (c m +i, Ch) there exist n — 1 linearly independent fields Vi(x, c) that 
generate a solvable Lie algebra and commute with v c (x). 

Let us add to Lagrangian (|7.3|) the term —x/2. Thus, we have placed the skate onto an inclined 
plane. Equations (|7.4|) - (|7.5|) hold if we replace x by h + x. Then the field Vh >c is equal to 

(\/ h — c 2 + x cosz, \/h — c 2 + x sinz, c). 

If h and c 7^ are fixed, then the fields 

vi = (2y/h - c 2 + x, -(cos z) /c, 0) , t& = (0, 1, 0) 

and Vfic are independent, and all their commutators vanish. In the same way one can solve the problem 
of rolling of a homogeneous disk on a rough plane, the problem of rolling of a ball in a vertical pipe 
and a series of other problems of nonholonomic mechanics. 



8. Existence of an invariant measure 

The existence of an integral invariant with a positive density is interesting not only from the standpoint 
of integration of differential equations. It is interesting by itself, from the standpoint of possible 
applications, for example, in ergodic theory. We are going to consider the problem of existence of an 
invariant measure for systems of differential equations. We are especially interested in its applications 
to nonholonomic mechanics. 

By the theorem on rectification of trajectories, in a sufficiently small neighborhood of an ordinary 
point there always exists an invariant measure with a smooth stationary density. Therefore, the 
problem of existence of an invariant measure is especially interesting near equilibriums as well as in 
sufficiently big domains of the phase space, where trajectories have the property of returning. Let 
us consider the first possibility. Let the point x = be an equilibrium of an analytical system of 
differential equations 

x = Ax + ... (8.1) 

We say that a set of (complex) eigenvalues Ai, . . . , A n of the matrix A is resonant, if m i^i = 
for some natural rrij. Note that a weaker resonance condition: ^mjAj = for some integer rrii ^ 
and Yl \ m i\ 7^ is usually used for investigation of system (|8.1|) (for example, in the theory of normal 
forms). 

Proposition 10. If a set Ai, . . . , A n is not resonant then in a small neighborhood of the point 
x = equation (|8.1j) has no integral invariant with an analytical density. 

The non-resonance condition is fulfilled, for example, in the case of ReAj ^ (^ 0) and 
£FieAi > (< 0). 

Proof. 

Let us expand the density M(x) in a convergent series with respect to homogeneous forms: 

M = M S + M s+1 + ..., s^0. 

Evidently, M s is the density of the integral invariant for the linear system x = Ax. One can assume 
that A as already reduced to the canonical Jordan form. Let us arrange the monomials of the form M s 
in some lexicographical order: 

M — \ a r mi r mn 

s ~ / j u mi...m„ J 'i ■■■• b n ■ 

m 1 + ... + m n — s 
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It is evident that divM s (Ax) is some form of the same degree. By equating its coefficients to zero, 
we get a linear homogeneous system of equations with respect to a mi .,. mn . The determinant of this 
system is equal to the product 

11 [(mi + l)Ai + . . . + (m„ + l)A n ] . 

7?l^ ^ 

m]_ + . . .-{-rrin — s 

This product is nonzero by supposition. Consequently, all a mi ...m n =0. ■ 

Remark. If a more strict condition of absence of resonant ratios in the traditional sense is fulfilled then 
equation l|8.1|l has no first integrals analytical in a neighborhood of the point x = 0. 

Let us consider as an example the problem on permanent rotations of a convex rigid body with 
an analytical convex bound on a horizontal absolutely rough plane (see [4]). The motion of such 
body is described by a system of six differential equations that have the integral of energy and the 
geometrical integral. In a particular case, when one of principal central axes of inertia of the body is 
orthogonal to its surface, we have a one-parameter family of stationary rotations about the vertical 
axis of inertia. Singular points of equations of motion correspond to the stationary motions. The 
characteristic equation has the following form: 

A 2 (a 4 A 4 + a 3 A 3 + a 2 A 2 + cqA 1 + a ) = 0. 



The dependence of the coefficients a s on numerous parameters of the problem is rather compli- 
cated; practically, they are arbitrary. The existence of the double zero root is connected with the 
existence of two independent integrals, since the differentials of the energy integral and of the geomet- 
rical integral are independent in the general case at points that correspond to permanent rotations. 
Fixing the levels of the first integrals we get differential equations on four-dimensional manifolds that 
have no invariant measure with an analytical density in the general case. Consequently, the initial 
equations also have no invariant measure in a neighborhood of stationary motions. 

Now let us consider the problem of existence of invariant measure for 
systems of differential equations that are similar to integrable systems that 
satisfy the conditions of theorem 1. It is natural to take the constants of the 
first integrals 1%, . . . , / n _2 as independent variables in a neighborhood of 
invariant tori of a non-disturbed integral system and to take angle variables 
x, y mod 2n on the invariant tori. In these variables, the perturbed system 
x has the following form: 




A(/) 



Fig. 3 



$(/, x, y) 



is = efs{I, x,y) + .. 
+ eX{I, x,y) + ..., 



1, 



n 



til) 



$(/, x, y) 



+ eY(I, x, y) 



(8.2) 

We assume that all functions in the right-hand sides of these differential 
equations are analytical in the direct product D xT 2 , where D is a domain 
in M. n ~ 2 = {Ii, . . . , I n -2}, T 2 = {x, y mod 2ir}; e is a small parameter. 
For system (|8.2|) . it is natural to consider the problem of existence of an invariant measure, the density 
of which is analytical with respect to /, x, y, 2-7r-periodic with respect to x, y and analytically depends 
on e: 



M = M + eMi + e 2 M 2 + ... 



(8.3) 



The unperturbed problem has an invariant measure with the density Mq. According to the 
well-known averaging principle, we average the right-hand sides of (|8.2|) with respect to the measure 
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dm = $ dx A dy. As a result, we get the closed system of equations for changing of slow variables / in 
the domain D 

i 8 =sF a (I), l<s<n-2, 

2tt 2tt 2tt 2tt 



(8.4) 





Proposition 11. Suppose mX(I) + nfj,(i) ^ in domain D for all integer m, n that are 
not equal to zero simultaneously. If averaged system (|8.4|) has no invariant measure with analytical 
density, then initial system H8.2|) also has no invariant measure with density (|8.3|) . 

System (|8.4|) is simpler than (|8.2|) : the sufficient condition of nonexistence of an invariant measure 
for H8.4|) is given by proposition 10. 

Proof of proposition 1 1 . 

Coefficients Mq and Mi of (|8.3j) satisfy equations 

dx $ dy $ 



E ^ + + l M ° Y = - + 4 f ) • < 86 > 

Since A//i is irrational for almost all I £ D, equation (|8.5|) implies Mq = r(/)$. Substituting this 
relation into (|8,6|) and averaging with respect to x, y we get the following equation: 

Y,w TFs = - (8 - 7) 

s s 

Consequently, T is the density of the integral invariant of ()8.4j) . It remains to show that T ^ 0. If it 
is not true then Mq = 0. But in this case the function Mi + EM2 + ■ ■ ■ is the density of an invariant 
measure for (|8.2|) If M\ = 0, this operation may be repeated once more. The proposition is proved. 

Remark. One can show that (under conditions of proposition 11) if averaged system (18. 4f) has no analytical 
first integral in D then initial system i|8.2|> has no integral that can be expressed as a series go + eg± + . . . with 
coefficients g s analytical in D x T 2 . 

Let us consider in more details the particular case, when n = 3. The index s may be omitted. 
Let F{I) ^ 0. If F{I) = at some point of the interval D then ()8.4|) evidently has no invariant 
measure. Therefore, we assume that F{I) 7^ in D. Consider the following Fourier expansions: 



F 



— = ^2X mn (I)expi(mx + ny), 



2^ Ymn(I) exp i(mx + ny), 

— = 2^ fmn{I) expi(mi + ny). 
The resonant set A is the set of points / £ D, such that 

E a mn 2 df mn . 1 -\t- \ 

— r— =00, a mn = — — + i(mX mn + nY mn ). 

I mX + na dl 

H + |n|^0 



Proposition 12. Suppose 

1) \(I)/n(T) const, 

2) the intersection Af) D is not empty. 

Then (|8.2|) /ias no integral invariant with density ()8.3j) . 
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Indeed, correlation (|8.7|) implies T = c/F, where c = const. Let 

= T,b mn (I) expi(mx + ny). 

Equation ()8.6|) gives us the set of correlations 

-(mX + n[i)b mn = ca mn . 

Let J € A. Then the condition 

£ \b mn \ 2 < oo 

implies c = 0. ■ 

The author is grateful to professor V. F. Zhuravlev who has read the paper and made a number 
of remarks. 
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